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In Gribov-Zwanziger scenario 1 - 1 of color confinement, Zwanziger proposed two types of 
horizon function in lattice Landau gauge, the first one 2 ' and the second one of which founda- 
tion was discussed in detail. 3 ' The second type horizon function is focussed on in the present 
study. Its derivation and the horizon condition are briefly reviewed along the line of Ref. 3), 
and it is also reviewed that this horizon condition and Kugo-Ojima color confinement crite- 
rion 6 ' coincide in the continuum limit. 8 '' 9 ' In case of free boundary condition in contrast to 
periodic boundary condition, it was shown that the horizon condition holds for each gauge 
fixed configuration in Landau gauge. 7 ' It is clarified that this fact can be derived from some 
identity relation which holds on an arbitrary network of links in Landau gauge. Thus the 
origin of the pointwise validity of horizon condition is highlighted. 

§1. Notations and generalities 

We intend to define Landau gauge SU{N) QCD on general networks, and to 
analyze standard lattice Landau gauge QCD from general point of view. So we 
have sites x and links £ = (x,x');pairs of sites in general networks. To each link I 
is assigned its intrinsic (positive) direction e M as convention, say e M = x' — x, for 
specification of basic link variable U e as a parallel transport U x y in the direction 
of e^, and of link- field (as "network current"), A x „ denote a value on the link 
£ = (x, x-\-eu) called as e M -component. Now links I with are considered as vectors 

xx'(= xx 1 ). We use notations U xx i = Ue = U x>e = U x ^ interchangeably, where 
x' — x = e^. Let l(x) denote a set of links, 

£(x) = {£\£ = xx =e il orl = x'x = e M } (1-1) 

and let £+(x) be defined set of positive link at x, 

e+(x) = {£\ i = xx' = e„} (1-2) 

and similarly £-{x) a set of negative link at x, 

£_{x) = {£\£ = x'x = e ll }. (1-3) 

It is to be noted that number of elements of £+(x) is not necessarily equal to that 
of £-{x) in general networks in contrast to periodic regular lattice or infinite regular 
lattice. 
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Thus gauge transformation by g G G is written for t G £+(x) as 

U$ = gJU t g x+eii , (1-4) 

and for t £ as 

Ujt = g x - e JU e g x . (1-5) 
We denote normalized antihermitian matrices A a as Lie algebra basis as 

[K,h]=fabAc (1-6) 

and 

(A a |A 6 )=tr(AtA 5 )=<5 ab . (1-7) 

All fields of the adjoint representation are often denoted as antihermitian fields in 
use of the above basis. 

We use bracket notation for suitable innerproducts for scalar fields (site function) 
and vector fields (link function), respectively, as 

<V#) = J>(^), (1-8) 

X 

(A^\B fl ) = J2^(A x jB x ^), (1-9) 

X,fJ, 

where we use a simple notation for link- functions, e.g., A^, and ^ implies summa- 
ry 

tion over x and £+(x), or over x and £-(x), which is equivalent to summation over all 
links y\ From site- fields (scalars) (j), two kinds of link-fields (vectors) are defined 

i 

such that 

d^<p = ^+e M - <f>x, (1T0) 
according to the associated positive direction of the link (x,x + e^), and 

^ = (4>x + e, + <f>x)/2, (Ml) 

where we use hereafter abbreviated notations x+fl for x + efJ . 

The following relations hold on any networks as well as on regular lattice. 

(A^) = (1-12) 
where the divergence of link-field is defined as 

{dnA^) x = ^2 A x ^ — ^2 ^x-n,n- (1"13) 

Adjoint of commutator reads as 

(A IX \[B IX C IX ]) = (-[B IX A„]\C„). (1-14) 
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Adjoint of *^ reads as 



where a site function is denned as 



^2 ^- x >v a x -^ 



> . 



(1-15) 



(1-16) 



§2. Definitions of gauge field and covariant derivative 

There are two possible options of A^(U), 
[/-linear definition; 

Ax,[m = (Ux,ii ~ f^r J ^)/2|tracclcss part) (2'1) 

log U definition; 4 ) 

U x ,„ = e A ^. (2-2) 

Next we define a covariant derivative D^QJ) which appears under an infinitesimal 
gauge transformation e £ as 5A Xjfi = D^(U)e] 

D^U)* = G(U^+[A^ (2-3) 



where the operation G(U^) on an antihermitian link variable, B^, is given by, 
[/-linear definition; 



1 \Up + Ul 



traceless part 



log U definition; 



with 



G{U„)B^ = S{A ll )B ll = {(^/2)/th(^/2)}S„ 



Ap = adj A „ = [Afj,,-]. 
It is to be noted that in both definitions G([/ M ) a 6 = G f (C/ /1 )b a from 



G(U,) ab = tv At- 



i \u u + m 



> h 



6a- 



and 



(A,)\ b = tr (At [A^, A b ]]) = tr ([A„ X a ][A„ A b ]]) = (^) 



6a' 



respectively. 

Adjoint of the covariant derivative is defined as 

(B^D^U)® = (-D^(U)B^\<P), 

where a site function D^([/)S M , with fi summation understood, is given as 



D„{U)B^ = d^GQJ^Bj + [A^B^f. 



(24) 

(2-5) 
(2-6) 

(2-7) 
(2-8) 

(2-9) 
(2-10) 
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§3. Optimization function and the Landau gauge 

The Landau gauge OA = can be characterized 5 ) such that 

SFu(g)=0 for v 5g, (3-1) 

in use of the optimizing functions Fjj(g) for each option of A^(U) as 
[/-linear definition; 

Fu(g) = £ tr {2 - (u^ + U°j) } , (3-2) 

x,n 

log U definition; 

Fu{g) = £tr (a« x J AQ ee {Al\Al). (3-3) 

x,n 

It is seen that in case of infinitesimal gauge transformations g~ 1 5g = e, variation 
5Fjj(g) is given in either definition as 

8F u {g)=2(A%\d ll e), (3-4) 

that is, by putting g = g(t), and g^g t = u t , 

^F u (g(t)) = 2(Af l \d,co t ). (3-5) 

It holds on the arbitrary networks that 

SFu(g) = -2(d^\s), (3-6) 

or 



^Fu(g(t)) = -2{d^\u t ), (3-7) 



d_ 

IE' 

which verifies the statement that a stationarity point of the optimizing function on 
a gauge orbit yields the Landau gauge. If we proceed further to a higher derivative 
in general, then we obtain that 

d2j? y )} = -2(8^(1/9)^) - 2{d^\uu). (3-8) 

Thus if Uf}^ = U^, then the variation of the optimizing function, AF = Fjj(e £ ) — 
Fu(l), is given up to the second order as 

AF = -2<VW + (e\ - d^(U)\e), (3-9) 

where it is put that e = rjAt with any site function rj constant with respect to t, i.e., 
uj tt = may be assumed in this case. 

Let us denote after Zwanziger, the Landau gauge space as 

r = {U\dA = 0}, (3-10) 
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and the Gribov region as 

n = {u\m{u) > o, u e ry, (3-n) 

where M(U) is a Faddeev- Popov operator, M = —dD(U), detailed properties of 
which will be investigated below. 

Now one can define a fundamental modular region A as a set of global 
minima of the optimizing function. 

A = {U\ Fu(l) < Fjj{g) for v 5 }. (3-12) 

It holds the following inclusion as 

a c n c r. (3-i3) 

So far, all notions and formula are valid on any networks, i.e., in any topology, 
and/or boundary conditions. 

For a while from now on, we assume d-dimensional regular lattice with L-periodic 
boundary conditions, and a set of link- field U on this lattice is denoted as L7 l, and 
a set of gauge transformation on this lattice, Gl- We define the iVth partial core 
of the fundamental region as 

A% = {U\ U G n L and ^(1) < F v (g), for v 5 € G NL }. (3-14) 

From Gl C Gnl, inclusion A^ C Al is easily understood, and one may write 

A% = A LN nn L , (3-15) 

and if is chosen to be a power of 2, say, then the partial cores are nested, 

A? CA% CA L for N' = 2 M ' >N = 2 M . (3-16) 
The core region is defined by the limiting set, 

E L = Af= lim Af = A OQ C\ n L . (3-17) 

In connection to the core region, the following theorem holds. 
Theorem 1 

Let U G Hi,-, where 77^ is a set of L-periodic configurations. The gauge transforma- 
tion g which brings U to Anl has a form such that 

g = he dx (3-18) 

where 

h G Gl, [^,^ = and e^ w = l. (3-19) 

Proof 

In the following, U is commonly used as denoting a configuration either being L- 
periodic or being N L-periodic. Let U 9 G A^l where g(x) G Gnl- Since the opti- 
mizing function F is an extensive quantity, it holds that the shifted configuration 
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in any direction by lattice unit have the same value of F. Particularly, the shifted 
configuration of U 9 in the negative ^-direction by L units is also a gauge transform 
of U into Anl, and it can differ from U 9 only by constant gauge transformaton, say, 
g^. Now since the shifted configuration U (x+Lfi) is identical to U(x) itself, it follows 
that g has a structure such that the L-unit shifted g, i.e., g(x + Lfi) in the negative 
^-direction is given by g{x)g il . It is obvious that [gfi,g v ] = 0, since translations in 
any directions commute each other, and gj^ = 1 from iVL-periodicity of g. We may 
write g^ = e s ^ L where O^s belong to the same Cartan subalgebra, and e s ^ NL = 1. 
Let h(x) = g(x)e~ ex , then it follows that 

h(x + Ln) = g(x + Lfj,)e- e » L e- 0x = g^g^e' 6 * = h(x), (3-20) 

and thus g = he 9x and h G Gl- q.e.d. 

Remark 2 
Let U G A%. It holds that 

Fu{l)<Fu{g) foT*g = he ex , (heG L ) (3-21) 

where 9^ = (M/NL)r/^ with the nonzero smallest elements rj^ of the Cartan subal- 
gebra such that = 1, and M is an arbitrary integer. It is to be noted here that 
U 9 G n NL , but U 9 (£ n L in general. 

Proof is selfevident from the definition of A^- 

Remark 3 
Let U G S L . It holds that 

Fu{l)<Fu{g) foT y g = he dx , (heG L ) (3-22) 

where 9^ = trj^ with the nonzero smallest elements rj^ of the Cartan subalgebra such 
that = 1, and t is an arbitrary real. Similarly to Remark 2, it is noted that 
U 9 G ilooi in general. 

Proof is selfevident from the definition of El. 

Now we investigate the behavior of the optimizing functions Fjj(g) under the 
gauge transformation relaxed so as to include the gauge transformation of the Bloch 
wave type. 

Let g be as g = he 0x = e^e dx = e^, and let us consider that 9^ and uj are some 
functions of t such that 

9^ = trjn and u> = to(t) (3-23) 

where r/^'s are suitably normalized constant elements of Cartan subalgebra, and uj{t) 
is a L-periodic scalar field with cj(0) = 0. Then putting 

h t = ^ = huj t and g t = ^=g^ t , (3-24) 
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we obtain 

£ t = e ~ ex Lo t e 9x + 7]x = cu' t + rjx, (3-25) 

where it is to be noted that neither ^ = ut nor ^ = £t hold in general. Another 
point to be emphasized here is that U 9 is considered as gauge transform of U G 
III C i^ooL an d is not of L-periodicity in general. It follows from general derivation 
before that 

d 



dt 

From (|3-25p . we have 



-F u (g(t)) = 2(Af l \d^ t ). (3-26) 



j t Fu(9(t)) = 2((A»\8^) + (A»\r,p)), (3-27) 

and 

j/v{g{t)) = 2(-(a /i i4«K> + (A^)). (3-28) 

The above equation (|3-28p is trivial as one on ooL-periodic lattice, but it has more 
implication than that. Although the periodicity of oj[ and A 9 ^ can not be demon- 
strated to be L, actual contribution from each link in (|3 27j) appears L-periodic, and 
the same holds in (|3-28p . and thus the derivation of (|3-28|) can be seen as such. 
Reasoning of this fact can be seen easily by noting that the inner product 
of link variables is invariant under the constant gauge transformation given 
by e ez+6,f1 / 2 i.e., the gauge transformation at the midpoint, and then there 
appear L-periodicity in the equation, and the subtraction can be inverted 
to the other side of the inner product with a minus sign. Explicit proof of 
this fact goes as follows. 

(AIM) = ((A h )f\d,(e- e *u t e e *)) (3-29) 

Let A X) n = A(U XjI j) denote the gauge field at a link (x, x + fi). Then under 
the gauge transformation of Bloch wave type, e 9x , it holds at each link 
that 

A e Z = = A(e'^ x+9 ^ 2 h +9 ^ 2 U x ^e +d -/ 2 e 9x+e ^ 2 ) (3-30) 

Then corresponding to the situation of constant gauge transformation, it 
reads that 

A ^ e -{ex+e^/2) e +e^/2 Ux e +e^/2 e ex+e^/2^ = e -(ex+e„/2) a(tj^ yex+e^/2 

(3-31) 

where \f x % = e +e ^ 2 U x ^e +e ^ 2 . And it holds that 
d^{e- dx u: t e ex ) = e -( fe +<W 2 ) 

(e-<VV+^ /2 " e e ^ 2 co t e~ e ^ 2 )e 9x+9 ^ 2 , (3-32) 
where ujt = 0Jt(x) and Wt,+u = ^t( x + A 4 )- Thus we obtain that 
((A h )f\d,(e- ex co t e ex )) = 

( e ^/ 2 (^e-<W 2 K+M> " (e- e ^ 2 {A h )^e e ^ 2 \u t ) (3-33) 
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where (A h ) x ^ = A((U h ) x ^) and it is to be noted that gauge fields ap- 
pearing in the inner proucts are L-periodic. Thus we can shift safely the 
expression as 



(3-34) 



where 



-11,11/ 



-ex 



{e e ^{A%^e- e ^) x ,, =A{e^Ut 

= e ex A(e-^- e ^Ut^e ex )e 
= e 6x A(U^)e-°*. 

Thus it holds that 

(e^\A h )^e-^\u t , + ,) = {e^A^e-^t) 

= (K,-M) 

Similarly it holds that 

(e-^iA^e^cot) = (Al\J t ) 
and thus we obtain that 

(Al\d,J t ) = -(d,Al\J t ) 
If we proceed further to a higher derivative in general, then we obtain as 

d2i y )} = 2(-(d,D,(U^M) - (d^Ku'At) + (D^U*)^)). (3-39) 

Now we consider a situation JJ 9 ^ G El and A 9 ^ = A^, and then since 

= D^U)(uo t + V x) = D„(U)uj t + G(U^ + [A x ^, VfM x% (3-40) 



(3-35) 

(3-36) 
(3-37) 
(3-38) 



±A9 



t=0 



the L-periodicity of A^ is easily violated due to the third term. Now we have 
d 



dt 



Fu{g{t)) 



t=o 



2(-(d^A lx \u t ) + (AM), 



(3-41) 



d 2 Fu(g(t)) 



dt 2 



t=o 



2(-(d li D lt (U)Zt\ut) ~ <VU(^)*} + {Dr(U)tt\ri»)) (3-42) 



where £t = ut + rjx, and then noting dA = 0, 



t=o 



d 2 Fu(g(t)) 



dt 2 



(3-43) 
(3-44) 



t=o 
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It should follow from the fact Un € yloo that 



Jt Fu( 9 (t)) 
d 2 Fu(g(t)) 



0. 



t=o 



dt 2 



> 0, 



t=o 



for v /i £ Gl and v r/^. For the estimate, 
we have 

(r/ M |LV([/)&) = (^I^CC/)^) + (r^G^)^) 

= -(^(C/H t K) + (r ?M |G(C/ M ) ??M }, 

where 

(t/JA^,^]} = (A x J[rJx^r]^}) = 
is used. From derivation, 



= d^D^{U)uo t + dpGifJ^ + [dpAp, r,x] + [3/> r? M ] 
= d^D^U)^ + d^G(U^ + 
= d lt D lt (U)u t + D li (U)ri li , 

we obtain for (|3-44p that 

1 d?Fu(g(t)) 



dt 2 



t=o 



(u t \ - d^D^{U)uj t ) - (utlD^Uftu) 
-(D^U) VlM \u; t ) + (^\G(U^), 



(3-45) 
(3-46) 

(3-47) 

(3-48) 
(3-49) 



(3-50) 



(3-51) 



1 d 2 Fu(g(t)) 

2 dt 2 



' D v \-dD\u t --^-D V ) 



t=o 



-3D 



-3D 



-(Drjl — lDrij + iri^GiU^) 



(3-52) 



Now we draw some necessary conditions for € from (|3-43p , (|3-45p , (|3-46p and 
([3^52]) . 

Theorem 2 

Let C/^ € El = n 77l. Then it follows that for the gauge transformation g = 
e wt e rjxt w j iere v w De i on g S to L-periodic scalar and ^t]^ belongs to the same Cartan 
subalgebra, the optimizing fucntion Fjj(g) behaves as 



d_ 

di 



Fu(g{t)) 



2(A M |»7 M ) = 0, 



(3-53) 



t=o 
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and 

1 d 2 Fu(g(t)) 1 1 

2 t=0 = ( " " " ^ " 

-(^IZ^I^) + (^IG(t^)M > 0- (3-54) 

Thus it is concluded that 

a^y, a ^ = ° ( 3 ' 55 ) 

X 

and 

(Drt\^\D v ) - {v^\G(U^) < 0. (3-56) 

Proof is selfevident from the fact that since Dr] is L-periodic for one can 
choose oj such that 

§4. Horizon function, horizon condition and Kugo-Ojima color 

confinement criterion 

Horizon tensor is defined as 

= -D^-dD^Dv - S^G(U^). (4-1) 

Taking the trace of the operator H^ u with respect to the normalized constant colored 
vectors ??^' a = 5^ u \ a with trA^A^ = 5 a b, one defines the horizon function H{U) as 

u,a 

= E« a l " D^-dD^DW) - (N 2 - l)E(U) 

u,a 

= h(U)V (4-2) 

where 

(N 2 - 1)E(U) = £ tv(X a tG(U Xtl ,)\ a ). (4-3) 

x,/i,a 

From Theorem 2, one has for U € that 

^ = ^-^^ = (4-4) 

X 

and 

#(£/) < 0, (4-5) 

where V = L A . 
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Zwanziger hypothesized that the dynamics on El tends to that on Al in the 
infinite volume limit, and derived the horizon condition, statistical average 

(h(U)) = 0, (4-6) 

in the infinite volume limit. Taking the Fourier transform of the tensor propagator 
of the color point source, 

(-D,(-dD)- l D v ) xayb , (4-7) 
one has, with assumption of the global color symmetry not broken, 

G^(p)5 ab = 5 ab [(e/d)(p^p u /p 2 ) - {V - {p^/p 2 )}u(p 2 )l (4-8) 

where 

e = (E(U))/V (4-9) 

and dimension d = 4. 

In local operator formalism of QCD, Kugo and Ojima proposed color confine- 
ment criterion 6 ^ based on the BRST(Becchi-Rouet-Stora-Tyutin) symmetry without 
Gribov's problem taken into account. Kugo-Ojima two-point function in the contin- 
uum theory is defined in the lattice Landau gauge QCD as 

(<V " ^KV) = ±Y, e ~ iP{X ~ y) (t^D,-^[A u ,X b ]) xy \ (4-10) 

where u ab (p 2 ) = 5 ab u{p 2 ) and it was shown that the sufficient condition of color 
confinement is given by 

u(0) = -1. (4-11) 

Putting Kugo-Ojima parameter as 

u(0) = -c (4-12) 

and comparing 

lirn G w (p) (4-13) 

with 

(h(U)) = (4-14) 
one finds that the horizon condition reduces to 

) = (£) +{d-l)c-e = (d-l)(c-?) = (d- l)h = 0. (4-15) 

Kugo-Ojima's and Zwanziger's arguments emerge to be consistent with each other 
provided the lattice covariant derivative naturally meets with the continuum one 

e/d = l. (4-16) 

This fact was pointed out and some numerical data were presented. 8 )' 9 ) 
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§5. Horizon condition in case of free boudary condition 

Let us consider the following quantity 

(d^lDrf) (5-1) 
on the arbitrary networks with assumptions that dA = in the sense of (jl-13p . i.e., 

(5-2) 



E w+- E 



where we put indices ± for link- field if /i € ^±(x), and with an assumption 
that a scalar function (p does not have the zero-eigenvalue eigenvector component of 
-dD(U). It follows that 

{d^Drf) = {4>\ - dD\4>) 



(5-3) 



One finds the Faddeev-Popov operator is symmetric when dA = which is seen 
below. 

D^<j> = d^G(U^d^ + [A^d^f, (54) 



(5-5) 



i-V'Vf -<uv"i)_ = E ([^^(^)]+-[^.^]+) 



/iG£_(x) 



1 



+ E M^ + ^Wl- 

/i££-(a!) / 



E ^ - E ^ 



Thus it hold that if d.4 = 0, then 



(5-6) 
(5-7) 
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and 

w zi^ \ Dd ® - M D ^) = °- ( 5 - 8 ) 

This equation (|5-8p is an identity which holds on any networks when dA = and (ft 
is free from zero-eigenvalue eigenvector component of —dD. 

Now let us consider the case with free boundary condition of regular lattice L d . 

and let us assume (ft = xr/ where r] p s are suitably normalized antihermitian matrices 
such that 

h,Vu} = 0. (5-9) 

Then it holds from 

d P <ft = ^ (5-10) 

and from 

{vM^]) = (MW,^]) = o (5-ii) 

that 

(Dm\^d\ D »V») - (Vp\G(U p )\rj p ) = 0. (5-12) 

By putting 

<° = S pp X a , (5-13) 

one obtains 

D^f = DuV P u' a = D p \ a (5-14) 
and then the vanishing horizon function 

H ( U ) = E ((D p \ a \-^\D p \ a ) - (\ a \G{U p )\\ a )\ = 0. (5-15) 

p,a ^ ' 

It is to be noted that A a in equation (|5-15p is located on links, and as seen from 
(|5-12p . D p acts on A a as defined in (|2-10p with non- vanishing first term. 



§6. Discussions and conclusions 



In §1, notations and generalities are given in order to discuss the case of free 
boundary condition from more general point of view, i.e., network QCD. The main 
purposes are to discuss behavior of the optimizing function for Landau gauge, gauge 
non-invariant function, and we need extended definitions, e.g., d p (ft and d p B p , with 
clear distinction between site-functions and link-functions, and finally to obtain an 
identity which holds on arbitrary networks. As a matter of course, full formulation 
of network QCD is out of scope of the present study. 

In §2, definitions of gauge field, [/-liner type and log [/-type, are given, together 
with covariant derivatives for each type, respectively, where difference between co- 
variant derivative (|2-3|) and covariant divergence (|2-1U|) should be noted. 

In §3, following Ref. 3), various kinds of regions in Landau gauge are defined, 
i.e., Gribov region, [2l, fundamental modular region, Al, and core region, El, on 
regular lattice of period L, where the following inclusions hold, 



e l = (/loo n n L ) c a l c n L . 



(6-1) 
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Theorem 2 states that for each J7„ € El = n ilz,, the horizon function defined 
in §4 (|4-2p takes non-positive value, H(U) < 0. 

In §4, although we have skipped the derivation of the horizon condition, (|4*6|) . 
Zwanziger showed in Ref. 3) that it can be derived from statistical average on aug- 
mented core region, !?£, in infinite volume limit, L — > oo, where 

"L C = £T(J7) < 0, C/ e tt L } c « r . (6-2) 

It is reviewed that the horizon condition and the Kugo-Ojima criterion of the color 
confinement 6 ) coincide with each other in the continuum limit. 8 )' 9 ) 

In §5, we focuss on the fact that the horizon condition holds for each configura- 
tion in Landau gauge on finite regular lattice with the free boundary condition. 7 ) It 
is found that the fact can be derived from the equation (|5-8j) . 

that is an identity which holds on any networks when dA = and <f> is free from zero- 
eigenvalue eigenvector component of —dD. Obvious reason why vanishing horizon 
function (|5T5p is not realized in the case of the periodic boundary condition is that 
(|5-10p does not hold. In a special case of U = 1 with free boundary condition, (|5T5p 
turns out to be 

H(U) = ((dpXa\^\d p X a ) - (A a |l|A a ) J = 0, (6-3) 

p,a ^ ' 

where there appear ±5 X 'charge density' only on the boundary surfaces. However, 
in the generic non-constant U cases, 'charge density' D p \ a in (15-151) spreads over 
d-dimensional volume. Thus one of the most important open questions is if the 
boundary condition affects the physics in the 4-dimensional bulk system in the ther- 
modynamic limit, i.e., continuum limit. 
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